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Abstract 

For an asymmetric double-well potential system, it is shown that, if the potential is quadratic 
until it reaches several times of the zero-point energies from the bottoms in each well, the energy 
eigenvalues of the low lying excited states of the double- well system must be close to the eigenvalues 
of the quadratic potentials. These eigenvalue structures suggest a method for the coherent control 
of the tunneling as well as realizing almost complete localization of the wave packet in one of 
the wells, by handling the double-well asymmetry. Numerical examples are included to indicate 
that the method could be useful also in a more general potential, and to propose experimental 
confirmations. 
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I. INTRODUCTION 



Quantum tunneling through a barrier is a fundamental physical effect , and the 
neutral atoms trapped in optical lattices have given us an opportunity to study the aspects of 
tunneling, including coherent dynamics on macroscopic scale. For the symmetric double-well 
system which has long served as a paradigm of quantum physics, the eigenvalue structures 
are well-known for the low lying excited states js-7|, and it has been pointed out that, 
by adding a specific driving force, the tunneling dynamics can be brought to a complete 
standstill known as coherent destruction of tunneling (CDT) {7]. Recently this coherent 
destruction has been visualized in single particle tunneling jsj as well as in the tunneling of 
Bose-Einstein condensates (BECs) [9|. On the other hand, in experiments, localized wave 
packets have been prepared in one well of an asymmetric double-well potential, and the 
tunneling dynamics has been observed by turning off the asymmetric part of the potential 
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Cs of interacting particles 



14J, to result in the non- 



are asymmetric in the asymmetric double-well potentials 
vanishing relative phase evolution rate [l5, 

In this article, we will find the energy eigenvalues of low lying excited states of an asym- 
metric double- well potential V(x) which has one local maximum at x = x c between the wells, 
by constructing WKB wave functions with the quadratic connection formula. To quantify 
the degree of tunneling, we define the tunneling visibility for a wave function ip(x,t), as 
V = [P m ax - Pmin) / (Pmax + Pmin) , where P max {P m i n ) denotes the maximum (minimum) 
value of P r {t) = ip*(x, t)ij)(x, t)dx during the time evolution. We find that the wave 
functions of (almost) arbitrary V ranging from to 1 can be realized from a Gaussian wave 
packet, by controlling the potential energy difference between the bottoms of the double 
well. The case of V ~ with P max ~ 1 amounts to the CDT, and this case can be realized 
when the two wells can be considered to be separate. A Gaussian wave packet at a stand 
still can be realized also in the well of higher bottom. 

While the results for a single particle can be applied only for the systems of noninteracting 
particles, we note that interacting systems have been extensively studied 17H20|. Particu- 



larly, for the systems of a few bosons, highly delayed pair tunneling analogous to nonlinear 



self-trapping has been found for the medium range of the interaction strength 



2l|. Though 



we only consider the double-well potentials bounded from below, the system of a particle 
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in a periodic potential with an additional constant force has been of great interest 
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24j, 
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and we note that this system has been analyzed through the instanton method 
is intimately related with the WKB analysis 

In the next section, before presenting the main analytic results, two asymmetric systems 
will be numerically studied. In Sec. Ill, we will construct the WKB wave functions for a 
general potential V(x). It will be shown that the asymmetric systems can be classified into 
two different regimes: In the one regime, an eigenfunction of low lying excited states has 
significant amplitude in both wells as in the symmetric systems, while, in the other, the 
eigenfunction describes the particle mostly localized in just one of the wells. In Sec. IV, 
we will develop formulas for the estimation of the energy eigenvalues in the regime of the 
localized eigenfunctions, and for the estimation of the tunneling visibility in the other. 
In Sec. V, the asymmetric double oscillator model will be numerically solved to give an 
implication on the eigenvalue structure of a general double-well system, and to show that 
WKB description could be remarkably accurate. The last section will be devoted to a 
summary and discussions. 



II. COHERENT CONTROL OF TUNNELING: NUMERICAL EXAMPLES 



In this section, we will study two systems numerically to indicate that the coherent control 
method could also be useful in a general potential whose wells are not exactly quadratic, 
and to expose that gravity may be used to control the tunneling dynamics in the settings of 
the recent experiments ja. \\\. 

First, we consider the system of a particle of mass m in the quartic double- well potential 
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with Iho = where C{a) is introduced to ensure that the minimum of the potential is 0. 

While the term proportional to a is added to give the asymmetry, Vq(0; x) is a special case 
of the well-known potentials js ?]. For Vq(0; x), the angular frequency for small oscillations 
at the bottoms of each well is u, and the barrier height is 3hoj/2. 

For the Gaussian wave packet </>G(a; x) =< x|</>g(oO >= exp[— (x — a(a)) 2 / (2lf)]/ (l^n) 1 / 4 
centered at the bottom of the right well, a(a), with /~ 4 = m = p ^^i^ \x=a(a), in 
Fig. [U we evaluate the probabilities Pi(a) = \ < <^G(a)|<^i(a) > I 2 ' where <pi(a;x) =< 



x\(fi(a) > are the eigenf unctions arranged in the order of ascending energy eigenvalue Ei(pt) 
(i — 0, 1, 2, • • ■). When we define the (unnormalized) wave function ?p a (x, t) as 



i/) a (x,t) 
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Fig. [T] indicates that, in this system of deep quantum regime, ip a {x, 0) closely describes the 
Gaussian wave function </>c(a; ^) for 0.01 < a in the given range of a. Further, plot (e) shows 
that, if 0.01 < a < 0.9, the Gaussian wave packet is closely described by an eigenfunction, 
ipi(a;x), which is in turn mostly localized in the right well of the higher bottom (We note 
that localized eigenfunctions have also been known in the buried double- well systems 261]) . 
As is well-known in the systems of symmetric potentials j^l, 5], the wave function 

, / -iE (0)t -iE 1 (Q)t\ 

( (po{0;x)e « + (fi(0;x)e n J describes the system which tunnels back and forth 
between an almost Gaussian state localized in the left well and the state localized in the 
right well, with V ~ 1. As, for 0.01 < a < 0.9, the ground state wave function (po(a; x) and 
the first excited state wave function <pi(a; x) are mostly localized in the left and right wells 
with the shapes close to a Gaussian, respectively, if we could change a from to a value 
which is larger than 0.01 but smaller than 0.9, without distorting the wave function by the 
change, then we will have a system where the probability density is almost stationary. For 
this stationary system, the probability of finding the particle in the left (or right) well de- 
pends on the details of changing a. If the change could be made in a much shorter period of 
time compared with the period of tunneling 2hir / (-E'i(O) — -E'o(O)), the probability of finding 
the particle in one of the wells crucially depends on the timing of the change. If we tune a 
back to or to 1, tunneling dynamics appears again, but in this time the tunneling visibility 
V could be less than 1 depending on the process. 

Second, for an atomic spinor trapped in a double-well of the potential 



V L (f3;x) = CE R [cos 2 kx + ^cos 2 (—)-^/2 + e/lQ} 



(3) 



+m(3gx, 

with E R = {hkf/2m flQ, we explicitly consider the case of C = 10 and £ = | which is 
similar to an experimental situation jg]. For (3 = 0, with the vanishing boundary condition 
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FIG. 1: (Color online) Tunneling dynamics in Vq{ol;x). (a) A potential (thick solid line) and 
three lowest eigenvalues (thin black, green, red solid lines), for a = ao = 0.985. The dotted 
and dashed, dashed, solid vertical lines indicate the values of a(ao), —b(ao), x c (cto), respectively, 
where Vq(«o! — b(ao)) = 0. (b) \(Pg{&o', x)\ 2 (solid line), and \ip ao (x, t)\ 2 at t = (dotted line) 
and at t = Jm/^E^ipto) — E\{q.q)) (dashed line), (c) and (d) The calculated tunneling visibility 
for i(; a (x,t) (solid lines) and an estimation (dotted lines). The estimation is made from Eq. (|16p . 
For (c), S e /5 a = 2V a (a(a)) / (Ei(0) — Eq(0)). For (d), since £2 (a) — E\{a) has the minimum at 
a = ai = 1.00, 2[V a (a(a)) - V ai (a(a 1 ))]/(E 2 (a 1 ) - ^i(ai)) is used as 6 e /5 a . The "+" mark 
denotes the value of ao- 
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at x = ^jr- (j :integer), there are three energy eigenvalues E\ (i = 0, 1, 2) under the barrier 
height, with E{ - Eft = 0.122^, E[ + Eq = 5.70E R , and E% = 7.27E R . 

If A = 2ir/k = 811 nm, g = 980 cm/s 2 and cesium atoms are in Vl(1]x), the fact 
mgX/2 = 0.580E R = 4.75(£'f / — Eq) then suggests that the tunneling visibility is very low 
for a cesium atom trapped in the vertical optical lattice aligned along the Earth's gravity. 
If the right minimum of a double- well is located at x r , we construct a Gaussian wave packet 
(/>(/; x) =< x\<j){f) >= exp[-(x - x r ) 2 /(2/ 2 / 2 )]/(/ 2 /» 1 / 4 , here, with a fitting factor /. 
Indeed, for the first excited state \(pf > of (3 = 1, we find that | < 0(0.759) \(p{ > | 2 = 0.980 
and | < 0(1.00) |</?f > | 2 = 0.945, which shows that the eigenfunction v?f(x) is closely 
described by a Gaussian wave packet, to prove the low visibility of the wave packet. As the 
visibility is high for the symmetric horizonal lattice, this shows that gravity may be used to 
control the tunneling dynamics in the optical lattices. 



III. WKB WAVE FUNCTIONS 



In this section, we will construct WKB wave functions for a general potential V(x), 

2 2 

assuming that V(x) is written as ^-(x + b) 2 and as ^^(x — a) 2 + ehuj r around the bottoms 
of the left and right wells, respectively. The eigenvalue structure will then be found by 
requiring the WKB wave functions to be asymptotically matched, in the overlapping regions, 
onto the exact solutions of the quadratic wells. 

In the quadratic regions of V(x), the eigenfunction is described by the parabolic cylinder 
function D v (z), and the eigenfunction of an energy eigenvalue hw r {y + e+ |)[= Hcui([i + |)] 
is written as 
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near the bottoms of the left and right wells, respectively, with U = \J {i = l, r )- On the 
other hand, by taking x c = 0, in the region of the barrier we have an approximate solution 
for the eigenfunction through the WKB method , as 
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FIG. 2: (Color online) An illustration of turning points. 



where p(y) = ^2m[V(y) — [y + e + ^)hu r }. The relations between the real coefficients 
Cl-i Cr, Nl, Nr may be given by comparing the eigenfunctions in the regions where the 
descriptions by the parabolic cylinder function and by the WKB function are both valid. 
On the negative real axis, the asymptotic expansion of the parabolic cylinder function is 
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The WKB wave function can also be expanded, for instance, as 



(6) 
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av p(y), 
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when ^> 1 and x is in the region of quadratic potential of the right well, where a u , —b. 



are turning points satisfying V(a u ) = V(—b^ 



e + \)fajj r (a > a v > 0, b > > 0) 



(see Fig. [2]). By comparing the leading terms in the parabolic cylinder function description 
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which originates from the first term in the right hand side (r.h.s.) of Eq. fl6]) with the relevant 
terms in the WKB description, we have 



r(-/i) V h \A t + \ 
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If any of fi or z/ is not close to a non-negative integer, in the large separation limit of 
e S\ P -ir d y > i and e / a " ^rfs/ > i ; Eq. © yields 

\N R \^>\C L \ and |iV L |>|C R |. (9) 

If the WKB condition p 2 (x) ^> is satisfied, in the region of the barrier, the first term 
in r.h.s. of Eq. (J5]) with positive (negative) Nr is a monotonically increasing (decreasing) 
function and the second term with positive (negative) iV^ is a monotonically decreasing 
(increasing) function. Eq. ()9]) then implies that if an eigenfunction exist for such u, it gives 
the probability distribution in which the probability of finding the particle in the barrier 
region is considerable. Since E v > J ip*(x)V(x)(p(x), if an eigenfunction gives a considerable 
probability in the barrier region, the eigenvalue can not be much smaller than V(0). 

On the other hand, if any of v or \x is close to a non-negative integer, the relations in 
Eq. (Q are not valid. If v is close to a non-negative integer, due to the singularity in the 
gamma function, the second term of the r.h.s. of Eq. (Q can also be a leading term. By 
comparing this type of leading term in the parabolic cylinder function description with the 
relevant term in the WKB wave function, we have a relation between Nr and Cr. By 
combining this relation with the first one in Eq. (jSJ), we have 
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If fi is not close to an integer, Eqs. (jSJ) and f fTUj) imply that the eigenfunction gives the 
probability distribution in which the particle is mostly found in the right well. 
If fi is close to a non-negative integer, we have the relation 



C L Ik V2ir 



Cr V K cos(/i7r)r(— v) \v+\ 



In this case, if v is not close to an integer, Eqs. flH]) and (11 ip imply that the eigenfunction 
gives the probability distribution of the particle mostly localized in the left well. 

For an eigenstate whose eigenvalue is much lower than the barrier height, the eigenvalue 
thus must be close to fouji(m+ |) or hu r (n + 1 + e) (n, m — 0, 1, 2, • ■ •), the eigenvalues of the 
quadratic potentials of the wells. Furthermore, the fact Dk(V2y) = e~ y2 ^ 2 Hk(y)/(V2) k for 
a non-negative integer k, implies that, if the eigenvalue is close to hw r (n + \ + e) (huJi(m + 
|)), the eigenfunction of the double-well system must be closely described by ^°(a; / r ; a;) 

b; If, x)) around the bottom of the right (left) well, with an eigenfunction of a simple 
harmonic oscillator ip^°(c; /; x) (= if fc (^) exp[— ]/ a/ y/nl2 k k\). 



IV. TWO DIFFERENT REGIMES 



The analysis of the previous section shows that an eigenfunction of the low lying excited 
states in the large separation limit has significant amplitude either in both wells or in just 
one of the wells. In this section, we will show that the eigenfunction of significant amplitude 
in both wells must be accompanied by another eigenfunction to form a doublet. As in the 
symmetric case, a linear combination of the doublet is responsible for the tunneling, and we 
will calculate the tunneling visibility for the combination. For the eigenfunctions localized 
in one of the wells, we will develop a formula for the energy eigenvalue estimation. 



A. Tunneling visibility 



For the case that both fi and v are close to integers m and n, respectively, we define 
5 e = e + n — ^-m, so that \5 e \ is equal or less than the minimum of \ and In this case, 
the corresponding eigenfunction gives considerable probabilities in both of the left and right 
wells, and /i and v should be written as fi = m + 5 M and v = n + 5 U with \6„\, \S e \ <ti 1. 



From the fact that 5 h 
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With 



when 5 tJ 



6 ± = -(-5 t ±^6f+^l), 
<5_, the eigenfunction is written as 

= (-l) n 5 a ^(-b; h;x) + Mgo; / r ; x) 

while the eigenfunction of S u = 5 + is 

-(-l)»M#(-&; Z,; x) + «°(a; Z r ; z) 



(13) 



(14) 



■0 (x) 



(15) 



The formal expression of 5± in terms of S £ and <5 a can be understood from the fact that, for 



<C 1, the tunneling dynamics is essentially described by that of a two-level system(TLS) 
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31] . If ipn°(a; l r ; x) is written as a linear combination of i/j (x) and ip + (x), the visibility 
of the linear combination is 



V = l/[1 



1A 
2 [ 5 n 



e \2l 



(16) 



For the visibility estimation in Fig. [Tfc) and (d), the parameter is determined from the 
consideration that, when 5 e = 0, the energy splitting is given as 2hw r 8 a . The fact that the 
tunneling dynamics significantly takes place for \8 e \ <C 1, even with m ^ n, may be closely 
related to the resonant enhancement of tunneling in the multiple- well structures 22|, |26j. As 
in the BEC loaded into an asymmetric double- well potential 12[ , if there are noninteracting 
N atoms in the ground state of V(x) of small 5 e (<^C 8 a ), the number difference between the 
left and right wells is proportional to S e . 



B. Energy eigenvalue estimation 

For an eigenfunction ip(x) of the system of V(x) with the energy eigenvalue E, we have 
the identity 

E - hw(n+ -) 

= go [V(x) - V r ho (x)] ^(x)jj h n °(a;l r ;x)dx 

where V r ho (x) = %p ^^f 1 . In numerical calculations, this identity may be efficiently used in 
estimating the energy eigenvalue of an eigenfunction which is close to ip^°(a; l r ;x). 
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As the visibility also implies, when 5 e is much larger than 5q, ip(x) of E ~ hw r (n + e + |) 
is mostly localized in the right well, and around the bottom it will be closely described by 
ipn°{a; l r ; x), to give ip^ pp (l r ; x) an approximation of ip(x) in this well. In the other regions, 
Eqs. (5,7) and the WKB method can be used to find ip^ pp (l r ; x). Eq. ( fTT|) then may be used 
to find a correction to u, as 

v — n + e 

This approximation of a localized eigenfunction (ALE) can also be made similarly, for the 
eigenstate of E « Hwi(m + |) which describes a probability distribution mostly localized in 
the left well. 



V. PRECISION TEST: ASYMMETRIC DOUBLE OSCILLATOR 



In application of the WKB method for a symmetric double-well potential, it is known 
that the energy splitting could be found accurately, if the (ground state) energy eigenvalue 
and thus the turning points are appropriately chosen Q|. If a well is quadratic with angular 
frequency u, then (j + ^)tuj (J: nonnegative integer) may be a good estimation for an energy 
eigenvalue. 

In order to check the accuracy of the formalism we have provided, avoiding the turning- 
point problem as much as possible, we consider the system of the asymmetric double oscil- 
lator potential 28] 

2 



V D {e;x) 



J hu 




{ V2l ho J 


Huj 





for x < 0, 
for x > 0. 



For this system, since both wells are exactly quadratic, the eigen 



(19) 



unctions are described 



by the parabolic cylinder functions on both sides of x = [28|, |29( and the continu- 
ities of the eigenfunction and its derivative at x = can be used to find the eigenval- 
ues EQ(a),Ei(a),---. As in Fig. |3l the calculations indeed show that, when a is a few 
times of Iho, the eigenvalues of the low lying excited states are close to tvjj{n + e + |) or 
Tujj{m + |). To expose that the estimation through the ALE is not valid when frw5 e is order 
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FIG. 3: (Color online) Energy eigenvalues for some lowest eigenstates of the system of Vd( £ ; %)■ e is 
chosen as 0.01 for (a) and (b), 0.3 for (c), and 0.5 for (d). Calculated values (Solid lines); estimated 
values through the ALEs (dotted lines) and through the approximation to a TLS (dashed lines). 
For the ALEs of the ground and second excited states, Eq. and the turning points a n , b n 
satisfying Vo(e; a n ) = Vo(e; — b n ) = (n + ^)hw are used with n = and n = 1, respectively; and for 
the first and third excited states, Eq. (fT0"|) and a n , b n of Vo^^n) = Vb(e; — b n ) = (n + e + ^)fria are 
used with n = and n = 1, respectively. In the approximation to a TLS, do and bo are determined 
from Vo(0;ao) = Vd(0;— bo) = \ftuj. 

of or smaller than the energy difference of the adjacent energy eigenstates, we add the ratio 
2ehu/(E^{a) - E^{a)) (dotted and dashed line) in Fig. E^b). 

When e is as large as 0.3, the ALE gives better results than the approximation to a 
TLS practically in the whole range where both methods are applicable [Fig. EJ^c)]. As the 
approximation to a TLS is suggested by the WKB method, Fig. [3] indeed shows that WKB 
description could be very accurate. Fig. [TJc) and (d) also suggests that this accuracy is 
not limited to the systems of the wells which are exactly quadratic. This with the reasons 
the WKB method provides implies that, if the potential V(x) is quadratic until it reaches 
several times of the zero-point energies hwi/2 and hu r /2 from the bottoms of the left and 
right well, respectively, the energy eigenvalues of the low lying excited states of the system 
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must be close to the eigenvalues of the quadratic potentials. 



VI. CONCLUSIONS AND OUTLOOK 

We have shown, through the WKB method of quadratic connection formula, that the sys- 
tems of asymmetric double-well potentials can be classified into two different regimes. In the 
regime of eigenf unctions giving significant amplitude in both wells, the tunneling dynamics 
could take place, while there is no tunneling in the regime of localized eigenf unctions. In 
this respect, the systems of the eigenfunctions mostly localized in just one of the wells are 
very different from those of the symmetric potentials. As Fig. [H(c), (d) and Fig. [3] clearly 
show, the WKB description could be very accurate, and the results given here may be valid 
for a system of the potential whose wells are not exactly quadratic. 

For the regime of localized eigenfunctions, even in the deep quantum limit, it may be pos- 
sible to confine a large number of noninteracting bosons in just one of the wells. For single- 
component fermions, in the light of the particle density pn(e; x, t) = Y\ ^^(e; x, t)ipja(e; x, t) 



(see, e.g., Ref. [32|), the number of fermions which can be confined in one of the wells is lim- 
ited by that of the localized eigenfunctions ipm(e] x, t). For the system of particles confined 
in just one of the wells, the tunneling dynamics can be initiated and controlled by adjusting 
e the potential energy difference between the bottoms of the double well, since, if we change 
e so that 5 e <C 1, ipmit', x, t) turns into a linear combination of the eigenfunctions of the new 
system. 

In the periodic arrangement of double-wells of the optical lattice where the tunneling is 
accompanied by a precession of the atom's angular momentum {J 10], a considerable time- 
periodic fluctuation of the population of atoms in a spin state could imply that the atomic 
spinors are in the states of high visibility. Since 5 a is very small in the large separation limit 
and the period of tunneling is inversely proportional to 5 a , if a tunneling phenomenon can 
be established over a long period of time, it can be used for precision measurements. 
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